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1 Introduction 



Quite often mass effects in high-energy collisions are considered as some not very spectacular 
corrections which finally die off. Nonetheless, it appears that in e + e~ annihilation even such 
overall characteristics as hadron multiplicities are quite sensitive to the value of masses of the 
primary qq pairs Q. 

Recent considerations have shown that calculations based on QCD agree well with the data 
at high enough energy [fj] and that they yield an asymptotically constant difference between 
multiplicities of hadrons induced by the primary quarks of different masses. 

In this paper we consider the possibility of a similar effect in a deeply inelastic process. 



2 Calculation of quark mass dependence 

Let us consider, for defmiteness, deep inelastic scattering of the electron (muon) off the proton. 
The hadronic tensor (an imaginary part of the virtual photon-proton amplitude) is defined via 
the electromagnetic current J M : 

W^(p,q) = \{^f J d A zexp(iqz){p\[J^z), J v {0)]\p), (1) 
where p is the momentum of the proton, p 2 = M 2 , and q is the momentum of virtual photon, 

g2 = _Q2 <Q 

A symmetric part of W^ v has two Lorentz structures: 

WW = (-9»„ + ^ F x {Q 2 ,x) + ±(p>~ <Z„P) {pu - F 2 (Q 2 ,x), (2) 

where the structure functions F% and F 2 depend on Q 2 and on the variable 

x = 91 (3) 

pq + ^(pq) 2 + Q 2 M 2 

In what follows we will analyse the structure function F 2 of deep inelastic scattering with 
open charm (beauty) production at small x. In this section we consider the case of one single 
quark loop with mass m q and electric charge e q . A general case will be discussed in Section 2. 

At small x a leading contribution to F 2 comes from one photon-gluon fusion subprocess Q : 

W»„ = J t^tj^jC'^ (?, k;m q )d aa > {k)dpp t {k)Y a '^ {k,p), (4) 

where k is the momentum of the virtual gluon, k 2 < 0. The tensor denotes an imaginary 
two gluon irreducible part of the photon-gluon amplitude, while T af3 describes a distribution 
of the gluon inside the proton. A quantity d a p is a tensor part of the gluonic propagator. 
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Let us choose an infinite momentum frame 

/ M 2 M 2 \ 

^={ P+ AP MP -AP)- (5) 

Then the gluon distribution T at3 has to be calculated in the axial gauge nA = with a gauge 
vector n M = (1, 0, 0, —1) [|J. One can take, for instance, 

*V = % + X P» ( 6 ) 

with x defined by Eq. (Q). 
From Eq. (0) we get 

Wi™ = F 2 {a) +F 2 (6) . (7) 



if, 



3g 



(pg) 2 + Q 2 M 2 



Two terms in the RHS of Eq. (0), F 2 and F 2 \ correspond to two tensor projectors, and 

Note that the structure function F L = F 2 — 2xF 1 is completely defined by the term p^p u 
and, thus, proportional to F 2 \ 

By definition, the gluon distribution T a/3 can be rewritten in the form 

r< * = r E S (P + fc - Pn)(p|/^(0)|n >< n|J|(0)|p>, (8) 

n 

where J* is the conserved current. Both \p) and |n) are on shell states that result in 

k a T a p = 0. (9) 

From an explicit form for (see Appendix I) one can verify that it obeys the same 
condition: 

fc-cgS = o. (io) 



Equations (Bj) and (110) allow us to simplify expression (13) and get (r = a,b): 



with the notations 



r<i a ) — — n n^ v 

,(6) _ __SQ 2 

(pq) 2 + Q 2 M 2 



C <*P ~ f^\2 , /02 AnP^P^aff- ( 12 ) 
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The tensor T al3 can be expanded in Lorentz structures 



k a kp\^ ( 1 pk\ ( 1 pk\ 1 

+ (*° ~ Ua L) - n ?L) A^ 3 + i Pa - Ha t) (** - Up ta) ¥ U (13) 

with I\ = Ti(k 2 } M 2 ,pk). 

Let us consider a contribution of the invariant function 1^ into the structure function F 2 (pT|). 
With the accounting for (|9|) and ( |10|) we obtain 

1 (r) _j)d* Q2 f X) df l-l 2 xyQ 2 z 2 (r) (Q^rr4x\ d a 

Z 2 9 i ^ 7 Z 2 1 + M 2 x 2 /Q 2 \ / 2 ' / 2 'J <91n/ 2 ^ ' J ' 1 } 



r a/3 = (^-^-)r 1 + (p a -A: a — JU,-^— J-r 2 



where 



and 



Here we used the notation: 



Z 2 = -k > 0, (15) 

z = - (16) 
pn 

M 2 z 2 

Ql = f^- (17) 



C (r) = -g^C^l 

In Eq. (0) the gluon distribution, g(l 2 ,z), is introduced: 

dr 2 



2(2 

If we use the new variable 



1 r dv 1 r 

g(i 2 ,z) = jjrj fkj.Wtk^z). (19) 

Ql 



—k 2 

t = 1 (20) 

pk + J{pk) 2 - k 2 M 2 



instead of kj_, we will arrive at the expression 



32tt 3 



2g 2 V 7 



sO 2 .*) = ^ / / ^ I M 2 + - I ri (/' 2 ,0- (21) 



A thorough analysis shows, however, that the main contribution to F 2 at small x comes 
from T 2 and T 4 in (|13|) (see also Appendix II). In Appendix I the following formula for F 2 is 
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obtained: 



r=a,b z 



2 (z—x)/x 



dp 
I 2 



+ Mr)(Q 2 m l A d ~ 2 



(22) 



As we are interested in a calculation of the difference of the structure functions correspond- 
ing to the massive and massless cases, we preserve those terms in which give a leading 
contribution to AF2. In Appendix I we have calculated the functions in lowest order in 
the strong coupling a s : 



C {a \u,v,y) 
C {a \u,v,y) 



where 



ct s 

Air 

ex. 

7Y 



{[(1 - y) 1 + y 2 ]L( Ul v, y) - [(1 - y) 2 + y 2 - 2v]M(v, y) - 1} 
y(l-y)M(v,y), 



L(u,v,y) = In 
M(v,y) -- 



u(l-y) 



y[v + y(l-y)] 

y(i - y) 



v + y{l -y) 



(23) 



(24) 



As for the gluon distributions, they are given by the formulae: 



G 



G 



32k/^/f«- Z »( M2 + F)t r ^« +r ^«l' 
Ql V 7 



r 2 i 
1 rdl' 2 



32<k 3 zJ Z' 4 



r 



v M2+ k 



(2g - z) 
4^2 



r 2 (r, + ^,0 



(25) 
(26) 



The analogous expressions for the functions are the following: 

C (b) (u,v,y) - 



= -^-y{2y[(l-2y)(l-y)-v]L(u,v,y) 

27T U 

o 

+ (1-y) [(1 - yf + y 2 - 2v]M(v, y)} + A/(l - y) 



2tt' 



12a, 1 



7T It 



y 2 (l-y) 2 M(«:2/)- 



(27) 



It may be shown that the leading contribution to AF2 comes from the region I 2 ~ m 2 , 
k 2 = —I 2 being the gluon virtuality (see Appendix I). Then one can easily see from fl24 ) and 
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(|28|) that the first two terms in C® are suppressed by the factor k 2 /Q 2 with respect to C^ a \ 
while the third terms in C® do not contribute to the difference C^\ m= Q — C^| m ^o- 

In the leading logarithmic approximation (LLA), only the function L remaines in Eqs. (p3"D , 
which results in 



(28) 



where P qg {z) is the Altarelli-Parisi splitting function and G(Q 2 ,z) is the gluon distribution in 
LLA defined by Eq. 



The gluon distribution fl25|) in LLA can be rewritten in the form (compare it with the 
corresponding formulae in Appendix II) 



(2 



(29) 



It is clear from (^) that AF 2 = F 2 \ m=0 — ^Im^o is defined by the quantities (r = a, b) 

AC^(u,v,y) = 0,y) - C^(u,v,y). 

By using Eq. (|23| ) we obtain the important result 

AC (a) = AC^\v,y) 
A(7 (a) = AC^{v,y), 



while from (|27|) we get 



AC® = -AC {b) {v,y), 



AC® = -AC®(v,y). 
u 



In this, we have 
So, we get || 



AC (a \ AC (a) \ 



-k 2 — >oo 



k 2 ' 



1 



X 



AF 2 {Q 2 ,m 2 ,x 



+oo e q 



1 <i,2 7 <i/ 2 



I 2 



AC 



m 2 x\ d 



+ AC 



I 2 ' z J d In I 2 
(m 2 n x\ d 



I 2 ' z d In/ 2 



G(Z 2 ,z) 
G(Z 2 ,z) 



Here 



at. 



AC(v,y) = ^\[(l-y) 2 + y 2 }ln 

47T 



1 + 



y(i -y) 



A<7(«,y) = ^(l-y)---^ 

7T u + y(l-y) 



(30) 



(31) 



(32) 
(33) 



(34) 



(35) 



5 



with G(P, z) and G(l 2 , z) being denned by Eqs. (|25|) and 



The integral in I 2 (|34"D converges because of condition fl3~3"f ). Contributions from AC^ and 
AC® are suppressed by the factors (m 2 /Q 2 ) InQ 2 and can thus be omitted. 



Let us consider the gluon distribution G At small z the leading contribution to G(l 2 , z) 
comes from the region z <C £, and we have 



G{l 2 ,z) ~G{l 2 ,z). 



(36) 



Taking expression (|36|) into account, the structure function F2 (^) has the following form at 
low x (with the term of the order of k 2 /Q 2 and m 2 /Q 2 subtracted) 



-F 2 = e 2 / - 



dz J(Z f X)/x dl 2 „ (Q 2 m\ x \ d 



l 2C Vl 2 ' l 2 > z) d\nl 2G{l2 ' Z) 



(37) 



where 



C{u, v, y) = - y) 2 + y 2 ]L( M , v,y)-[{\- 3y) 2 - 3y 2 - 2v}M(v, y) - 1}. (38) 



As for the difference of the structure function, we obtain the following prediction 



1 . „ ,■<, o s 1 . „ . n , o } dz J dl 2 K ^ ( m 2 x\ d 
-AF 2 (Q 2 ,m 2 q ,x) = -AF 2 (m 2 q ,x) = e\ j - j -^AC Uf, - 

Ql V ' 



<91n/ ; 



: G(Z 2 ,z), (39) 



where 



0:.s 



AC<(^) = ^[(l-y) 2 + ^ln 



1 + 



1/(1-2/) 



-(l-2y) s 



v + y(l -y) 



(40) 



3 Relation between measurable structure functions 

Up to now, we considered those contributions to F 2 that came from the quark with electric 
charge e q and mass m q , F 2 \ m ^Q. Then we have taken the analogous contributions from the 
massless quark with the same e q , -F2|m=o ; an d have calculated the quantity -F2|m=o — F 2 \ m ^ . 

The total structure function F 2 has the form 

F 2 (Q 2 ,x)=Y,e 2 q F«(Q 2 ,x), (41) 

where the functions F% are introduced (q = u, d, s, c, b). 
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The structure functions describing open charm and bottom production in DIS, F 2 and F 2 
respectively, are related to F% and F 2 by the formulae 

4 ~ 

2 g 2; 

F\ = \F b 2 . (42) 

At low x one can put (m u = = m s = is assumed) 

F% = F* = F° = F 2 (43) 
and define the difference between heavy and light flavour contributions to F 2 : 

AF 2 C = F 2 - F 2 C , 

AF 2 fe = F 2 - F 2 b . (44) 

Notice that there are the functions F 2 and F 2 that have been calculated in the previous section 
(see Eqs (|37|) and (g)). 

Let us now represent the function F 2 ([37]) in the following form 

1 Y 







where we denote 



2 



(46) 



and introduce the variable rj = ln(k 2 /QQ) 



Analogously, we get from 

x J y J alnr 



1 AifJ= fd V AC( Vl y)-J- 1 ,G(Y m -ri 1 -) 1 (47) 



X — oo 



with 
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F m = ln-§. (48) 

Here ?7 = ln(mg/fc 2 ?/(l — y)) ~ \n{m? q /k 2 y) (remember that we consider small x). 

The expression for AC is given by Eq. ([H]) and, in terms of the variables rj and y, looks like 



- 2 , 2n 



ln(l + e")-(l-2?/) 2 



1 + 



(49) 
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As for the expression for C, it has to be defined via relation ([TJ]) and exact formulae QL14J 
and ( |1.21| ) taken at m — 0. The result of our calculations is of the form 



1 . 1 + Z7 



— V + V 

u 2 



where 



—V 

u 2 



(50) 



u 

V 



l-4y(l-y)e-<i, 
(1-y) [y + (l-y)e-*] (l - e"") . 



(51) 



It is clear from ( |49|) that 

AC( V ,y)>0 (52) 
for — oo < < oo, < y < 1 and AC(rj,y) is negligible at r] < (see Figs. la-Id). 

Moreover, the quantitative analysis shows that at most at y < 0.2, which is relevant for 
small under consideration, one has 

C( V ,y)>AC( V ,y), V > 0, (53) 

(see Figs. 2a-2d). Neglecting the small contribution to F 2 from the region rj < and taking 
into account that dG(Q 2 , x)/d In Q 2 > at small x (cf. 0), we thus conclude 

AF$(m 2 q ,x) < F 2 (Q 2 ,x)\ Q 2 =m2q . (54) 



In terms of the observables F 2 , F 2 and F 2 , the inequalities (j5^) and ( p4|) can be cast in the 
forms 

(F 2 -2.5F 2 c -F 2 fe )(Q 2 ,x)| largcQ2 > 0, 
(F 2 -F 2 c -7F 2 6 )(Q 2 ,x)| largeQ2 > (55) 

and 

(F 2 -2.5F 2 C -F|)(g 2 ,a;) < (f 2 - F 2 C - F 2 b ) (Q 2 , x)| Q2=m 2, 
(F 2 -F 2 c -7F 2 b )(g 2 ,x) < (F 2 -F 2 c -F 2 b )(g 2 ,x)| Q2=m? . (56) 



Data on the total structure function F 2 for Q 2 between 1.5 GeV 2 and 5000 GeV 2 and x 
between 3 x 10~ 5 and 0.32 are now available 0. As for the charm structure function, there 
are preliminary low-x data on F 2 at Q 2 = 13 GeV 2 , 23 GeV 2 and 50 GeV 2 with rather large 
errors |T7fl. 



Using the first of the inequalities ( [56"D we get (assuming F 2 (m 2 
F 2 c (Q 2 ,x) > 0.4 \F 2 (Q 2 ,x) - F 2 6 (Q 2 ,x) 



c ,x), F 2 



[ml, x) 



F 2 (m 2 c ,x; 



(cf. Q)) 
(57) 
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Let us estimate F^Q 2 , x) from below for x = 5 x 10 3 and x = 5 x 10 4 and several values 
of Q 2 . There are data on F 2 (Q 2 ,x) for Q 2 = 2.5 GeV 2 , x = 4x 10~ 3 and Q 2 = 2.5 GeV 2 , 
x = 6.3 x 10~ 4 @. The relative bottom contribution, F\jF 2 ^ reaches at most 2 to 3% at HERA. 

Putting F 2 (m 2 , 4.0 x 10" 3 ) ~ F 2 (m 2 , 5.0 x 10~ 3 ) and F 2 (m 2 , 6.3 x 10~ 4 ) ~ F 2 (m 2 , 5.0 x 10" 4 ) 
and choosing m c = 1.58 GeV we obtain from (??) the lower bounds for F% presented in Tables 1 
and 2. 



Q 2 , GeV 2 


12 


25 


45 


FZ(Q 2 ,x) 


0.106 ±0.054 


0.139 ±0.044 


0.195 ±0.046 



Table 1: The lower bounds on F 2 c (Q 2 ,x) for x = 5 x 10 3 . 



Q 2 , GeV 2 


12 


25 


Fi(Q 2 ,x) 


0.207 ±0.063 


0.355 ±0.070 



Table 2: The lower bounds on F£(Q 2 , x) for x = 5 x 10 4 . 

These quantitative estimates of _F 2 C do not contradict the preliminary data on the charm con- 
tribution to F 2 0. For a detailed comparison of our predictions with the data, an improved 
measurement of the charm component _F 2 C is required. 

Conclusions 

In this paper we have demonstrated that the lowest-order quark loop contributions to the 
structure functions at small x contain mass-dependent terms which scale at high Q 2 . This 
effect can be observed experimentally, and we predict theoretical bounds for the corresponding 
contributions from c- and 6-quarks (see Eqs. ( |55|) and fl56"D , Tables 1 and 2). 
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Appendix I 



Ilp= - 9a/3-[(l 4 -2l 2 qk-k 2 Q 2 ) + 2m 2 k 2 }+r a r l3 2(Q 2 -2m 2 ) 



In this section we calculate an imaginary part of the photon-gluon amplitude averaged in the 
photon Lorentz indices with the tensors g^ v and p^p u (see Eqs. (|I2D). 

Let us consider the function C^}. In the first order in the strong coupling a s (a one-loop 
approximation), it can be represented in the form 

CapiQ 2 , k<2 , qk) = ^ J d 4 r[^I^{r, q, k; m 2 ) + — — I— — J^( r , q, k, ; m 2 )} 

x 5 + ((g-r) 2 -m 2 )<5+((A; + r) 2 -m 2 ), (1.1) 

where 

+ {r a kp + k a rp)(Q 2 - 2m 2 ) - (r a qp + r^)/ 2 

- + k a q^)l 2 (1.2) 

corresponds to a contribution of a ladder diagram, while 

lip = g a p\k 2 [{k 2 -Q 2 + 2qk)-2m 2 ]+r a r p 2[{Q 2 + k 2 ) - 2m 2 ] 

- k a k(i{l 2 - Q 2 ) + q a qp(l 2 + k 2 ) 

- (r a k p + Kr p ) [(I 2 -Q 2 -k 2 - qk) + 2m 2 ] 

- {r a q p + r a q p ) [{I 2 + Q 2 - k 2 - qk) - 2m 2 ] 

- {k a qp + k a q^[{Q 2 + k 2 )-2m 2 ] (1.3) 



is a contribution of a crossed one. In Eqs. (|L2|) and (]L3|) we denote 

Z 2 = m 2 - r 2 . (1.4) 

The following equality thus takes place 

d 4 rr a f(r 2 )5 + ((q — r) 2 — m 2 )5 + ((k + r) 2 — m 2 ) 

:[/c Q A!(Q 2 , /c 2 , gfc; m 2 ) + q a A 2 {Q 2 , k 2 , qk; m 2 )], (1.5) 



where 



A 1 = --Ljdl 2 [Q 2 (qk + k 2 )+l 2 (qk-Q 2 )]f(^ 



l 2 _ 
l 2 



M = 2^1 ' dl 2 [l 2 {qk + k 2 )-k 2 {qk-Q 2 )]f{l 2 ) (1-6) 
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with 



72 
f ± 

D 

s 



(5*0 ± 



£>(1 



4m 2 



1/2 



(g£;) 2 + Q 2 £; 2 , 
A; 2 - g 2 + 2gA;. 



Analogously, we have 



where 



d rr a rpf(r )S+((q — r) — m )5+((k + r) — m ) 

7T 1 

—=^5! + k a k B 2 + g^^s + -(A^ + q a k p )B 4 ], 
4v iJ ^ 



Bi = J dl2 {gD Rs+ 2 m2 } f{l2) > 
i 2 



Bo = dl 



Bo 



l 2 



dl 2 {^Rs + ±[(l 2 + k 2 ) 2 + 2m 2 k 2 ]]f(l 2 ), 



B A 



dl 



i 2 



\-W* qkRa + ^D [R ~ l2s ~ 2m2 ^} /(/2) ' 



In Eqs. ( p9|) a notation 

R = l 4 - 2l 2 qk - k 2 Q 2 

is introduced. 

Accounting for all that was said above, we get 



C$(Q 2 ,k 2 ,qk) 



dl 2 



^«/?(^ 2 5 5) k\ m 2 ) + j^-I^^l 2 , 5) k, 



I 4 



P(2qk-l 2 Y 



with 



1 af3 



1 - 



Q 2 

R — k a k/3^r 



D 



R-Q 2 s- —Q 2 sR 



+ q a qp 



D 



Rk 2 + l 4 s - ^k 2 Q 2 sR 
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and 



- (k a q/3 + q a kp) — 



Rqk - l 2 Q 2 s - —Q 2 qksR 



+ k a kp^ 
1 

+ QaQPp 



(k 2 - Q 2 ) + ^sR 
R-Q 2 s + Q 4 - ^Q 2 sR 



( /2 + fc 2 )2 _ k 2 Q 2 + ^_ k 2 gR 



D 



{Q z k z qk 1 

g a p - k a k p — + q a qp— - (k a qp + q a kp)—> 



R-{1 2 + k 2 ){qk - Q 2 ) + Q 2 qk - ^qksR 
k 2 



qk 



(1.12) 



f 2 



9a/3 

+ k a k g — 



1 



k 2 + ^(Q 2 + k 2 )R 



R-Q 2 (Q 2 + k 2 )-^Q 2 (Q 2 + k 2 )R 
R-k 2 (Q 2 + k 2 )-^k 2 (Q 2 + k 2 )R 



(k a q/3 + q a kp)—qk(Q 2 + k 2 ) 



1 + R 

2D 



- rn \ -g a/3 
1 



Q 2 - —sR 
v 2D 



k a k/3 



D 



R-2Q 2 {qk-Q 2 )-—Q 2 sR 



+ q a qp 



D 



R + 2k 2 (qk-Q 2 ) + —k 2 sR 



+ (kaq/3 + q a kp) 



D 



R - (qk) 2 + 2qkQ 2 + Q 2 k 2 - T^qksR 



( Q 2 k 2 qk 

2rn 4 lg a/3 - k a kp— + q a qp— - {k a qp + q a kp) — 



(1.13) 



Equation (|I.11| ) can be represented in the form 



f 2 



a 



a/3 



4:71 y/D 



p 



+ 



1 

p 
1 

2P 
2 
P 



Q 2 



2D 
1 



l - — Q 2 s 
2D^ 



k 2 Q 2 - 



P 



^- k (k 2 -Q 2 )(2(qk) 2 + k 2 Q 2 )s 




m 2 2 m 4 




2 > 


qk P qk 




'k 2 


^ 2D ^ 


m 2 



1 



1 



m 



2D 



k 2 s 



12 



i 2 



+ B, 



+ 
+ 



a/3 



a s 1 
4vr v^D 



/- 



1 

P 
1 

2^ 
2 



3 



- M + 3Q 2 + ^(2(gA;) 2 + fc 2 Q 2 )s 



m 2 m 
qk I 2 qk 



1 - w** 



k 2 Q 2 - 



l 4 



k 2 -Q 2 - ^k 2 Q 2 s 



m 



m 



1 + 



(1.14) 



where 



5, 



9 a/3 



k a kj3 
qk, 



qk. k 2 



(1.15) 



From Eqs. (p.!4| ) and (|I.15| ) one can see that obeys a condition 



fc°Cg = 0. 



(1.16) 



In deriving (|Ll4|) from (pI|)- (p3D we took into account that l 2 ± -> l\ at I 2 -> 2qk - I 2 , l 2 ± 
being defined by Eq. (|LT 



Note that the difference of the structure function under consideration, AF 2 , is defined by 
AC^g (r = a, b). The analysis shows that a leading contribution into AC^ comes from the 
region —k 2 ~ m 2 <C Q 2 . Rejecting terms 0(k 2 /Q 2 ) and 0(m 2 /Q 2 ), we obtain the formula 



(a) 
a/3 



J2 
Z 2 



1 



1 



1 



I 2 V 2(qk) 2 



Q 2 s) + 



1 fc 2 Q 2 
2/ 1 



1 - 



2(gfc)2 g2s + 2 F 



+ 



2(<7*0 

i fc 2 g 2 

2F qk 



+ B a 



a 



An 



s 1 dl 2 



ll 
3 



Ifl-^-^ 

Z 2 V 2(g/c) 2 



2(?fc) 



7 g2 s + 2 — 



1 



fc 2 / 2(gfc) 



(1.17) 



Now let us consider the function C^l (see Eqs. (0)). The result of our calculations is of 
the form 



C { 3(Q 2 ,k 2 ,qk) = ^/^[ljyr,g,fc;m 2 ) + 



2tt 2 D J ~ ' |/ ; ' ' ' l 2 (2qk-l 2 ) Jaf3i - r,q,k, ' ,m ^ 

x 5 + ((<? - r) 2 - m 2 )5 + ((k + r) 2 - m 2 ), 



(1.18) 
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where the contributions from ladder and crossed diagrams are, respectively, 

1 



g a p~k A s + Ta rp2[R - k 2 (2l 2 - Q 2 ) + k 2 s)\ 

+ k a kfi[R-k 2 {l 2 -g 2 )] 

+ (r a kp + k a r p )[2R - k 2 (3l 2 - 2Q 2 ) + k 2 s] 

+ {r a qp + r a qp)l 2 k 2 + (k a gp + k a qp)l 2 h 2 



(1.19) 



and 



J. 



g a p-k A s + r a r p 2(R - k A ) + k a k p {l 2 - Q 2 )k 2 

- q a qp(l 2 + k 2 )k 2 + (r a k p + k a rp)[R + k 2 (l 2 - k 2 - qk)} 

- (r a qp + r a qp)[R-k 2 (l 2 + k 2 - qk)} 

- (k a q p + k a q^[2R-k 2 (k 2 -Q 2 )}. 



(1.20) 



Using formulae ( p5|) and (|L8|) we obtain 



Mb) 



+ 
+ 



+ 



3a, Q 2 



1 

D 
1 

Z 2 
1 



l 2 _ 

,2 on2 



qk(3k 2 - 3Q 2 + 4qk) - 2k 2 Q 2 

2(k 2 + qk) + \k 2 (k 2 + Q 2 
qk 



k 2 s 



m 



l - — Q 2 s 
2D^ 



k 4 s — —k 2 sm 2 
Z 4 



+ l 2 ^(k 2 + qk)s-^ 3 (k 2 + qk)(k 2 + 2qk)s 



3a s Q 2 

a(3 87T ,03/2 



I 2 



li 

,2 on2 



-Af - Q 2 ) 

qk 



3 

D 
1 

P 
1 

2Z 4 



qk(3k 2 - 3Q 2 + 4qk) - 2k 2 Q 2 

2{k 2 + qk) + \k 2 {k 2 + Q 2 ) ni- 
gh 



k 2 s 



1 - ^ 



(1.21) 
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At -k 2 ~ m 2 < Q 2 we get from (|I.21|) 



a 



i 2 

( b) 4 3a s Q 2 f 2 fl 



(g/c — Q 2 )s 2m 



(gA;) 1 



+ 



A; 2 



1 A; 4 s 
2/ 1 



Q 2 s 2m 2 



3a 5 Q 2 
Q/3 8vr (gA;) 2 



2(g£;) 2 



k 2 



S i2 S 

h I 

qk 2(qk)' 



dV 



I 2 



I 2 



(3qk - 2Q 2 )s 2m 



l_k*s_ 
W^k 



3Q 2 s 2m 2 



2(qky 



k 2 



(qkf 



2s 



+ 



k 2 



3s 



qk + ^ 2(qky 



(1.22) 



When calculating C^l we used the replacement 



krv k 



(1.23) 

'{pq) 2 + Q 2 M 2 ^J(qk) 2 + Q 2 k 2 
This results in power corrections that are, however, ignored everywhere in our consideration. 
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Appendix II 



In anology with quark distribution, scalar gluon distribution inside the nucleon, D g , can be 
defined by considering deep inelastic scattering with scalar gauge-invariant gluonic currents 



{111) 



Let us define D g via a gluonic structure function F: 

D(Q 2 ,x) = -F(Q 2 ,x) = -3-DiscT, 
x Air 



where 



T = iJ d 4 zexp(iqz)(p\TJ(z)J(0)\p). 



Function F is given by the formula 



x An J (2tt^ 4 



d ' k l U^(q,k)^(k,p) 



Here U a p is a gluonic partonometer 

n Q/3 = -g^2n5((q + k) 2 )Q 2 , 



(11.2) 
(11.3) 

(II.4) 
(II.5) 



while T al3 means an imaginary part of the virtual gluon-nucleon amplitude with tensor parts 
of the gluonic propagators d aa > included. So, we have by definition 



T«p = d Q *(k)dpp(k)T o/f * 
where the gluon distribution enters Eq. (|). 
Tensor T a p has the following Lorentz structure 



(II.6) 



a/3 



d aP (p)Y l + d af3 (k)T 2 + [Pa-k, 



pn 

1 kn 



Pf3 ~ K/3 



pn\ 1 ~ 
ka-r- T^r 3 + 



kn) k 2 



3 + n a np 



(kn 



(n.7) 



with Tj = Ti(k 2 , M 2 ,pk). Then from Eqs. ( [II. 6|) , ( P. 7] ) and (13) the relation between Tj and Tj 
can easily be obtained. In particular, we have 



r 2 + r 4 



(II.8) 



If we substitute flTTSD , (|TTT7|) into (0X4]), we get in LLA 

? dk 2 r d 2 k 



D(Q 2 ,Qix) 



a A; 4 J 2(2n) 



T^k 2 , k±, x) + T 2 (k 2 , k ± , x) + -^^ T 3,(k 2 , k±, x) 



. (II.9) 



As can be seen, at small x the gluon distribution ( |11.9| ) is mainly defined by the invariant 
function 1^ (or, equivalently, by the combination T 2 + I^). 



16 



References 



[1] J. Chrin, Proc. of the 27th International Conference on High Energy Physics, Glasgow, 1994 
(Eds. P.J. Bussey and I.G. Knowles), p. 893; A. De Angelis, Talk at the EPS Conference 
on High Energy Physics, Brussels, 1995 (to appear in the Proceedings); W.J. Metzger, ibid. 

[2] V.A. Petrov and A.V. Kisselev, Z. Phys. C66 (1995) 453; Nucl. Phys. B (Proc. Suppl.) 39B, 
C (1995) 364. Early qualitative indications on the energy independence of the multiplicity 
difference were made by Yu. Dokshitzer, V.A. Khoze and S.I. Troyan, J. Phys. G 17 (1991) 
1481, 1602. 

[3] A.V. Kisselev and V.A. Petrov, Preprint IHEP 95-115, Protvino, 1995; Talk at the EPS 
Conference on High Energy Physics, Brussels, 1995 (to appear in the Proceedings). 

[4] S. Catani, M. Ciafaloni and F. Hautmann, Nucl. Phys. B366 (1991) 135. 

[5] J. Kwiecihski, A.D. Martin and P.J. Sutton, Preprint DTP/96/02, Durham, 1996. 

[6] HI Collaboration, S. Aid et al., Preprint DESY 96-039, Hamburg, 1996. 

[7] A. De Roeck, Talk at the 2nd Krakow Epiphany Conference, Krakow, 1996 (to appear in 
the Proceedings). 

[8] A.D. Martin, R.G. Roberts and W.J. Stirling, Phys. Rev. D50 (1994) 6734. 



17 



Figure Captions 

Figs. la-Id: AC(r),y) as a function of the variable rj at several fixed values of y. 

Figs. 2a-2d: C(r),y) (continuous curves) and AC(rj,y) (dashed curves) as functions of the 
variable r\ (rj > 0) at several fixed values of y. 
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